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ABSTRACT 



The use of the dual graph in determining the maximal flow capacity 
of an undirected source-sink planar network has been extended to general 
networks. A directed dual graph is first defined for directed source- 
sink planar networks such that the length of the shortest route through 
this dual is equal to the maximal flow capacity of the primal. The 
construction of the dual is then modified to handle flows in networks 
which are not source-sink planar. An algorithm is presented for finding 
the modified shortest route through this dual having a length equal to 
the maximal flow capacity of the primal network. 
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1. INTRODUCTION 



Ford and Fulkerson (1956) observed that an easy way to find the 

\ 

value of the maximal possible flow through an undirected source-sink 
* 

planar network is to construct the dual graph, assign the capacities 
of the intersected primal arcs as lengths of the undirected arcs in 
the dual, and then find the shortest route through the dual. The 
length of this shortest route will be equal to the value of the mini- 
mum cut set of the primal network and hence is equal to the value of 
the maximal flow because of the min-cut max-flow theorem. 

Extension of this idea to directed networks began with de 
Ghellinck (1961). He developed a convention for constructing a dual 
of a directed source-sink planar network but was not concerned with 
arc lengths. In a paper written in 1969, Doulliez and Rao (1971) 
provide the dual construction and the arc length conventions for 

i 

directed capacitated flow networks when zero lower bounds exist on 
arc flows. This convention was apparently the result of conversations 
with de Ghellinck. 

Because of a need for dual graphs of directed network in studies 
of network interdiction models, the author (1970) independently devel- 
oped the same convention as Doulliez and Rao but went one step further 
and provided a convention for handling problems where the lower bound 
on arc flow could take any negative value or any positive value not 
exceeding the upper bound capacity. This paper also observed that for 

Ford and Fulkerson (1956) defined source-sink planar networks as planar 
networks which remain planar after an arc connecting the source and sink 
is added to the network. 
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flow networks having positive lower bounds on arc flows, an infeasible 
primal problem is indicated by a cycle of negative length in the dual 
network. In a later paper (1971) the author proved that a cycle of 
negative length was a necessary and sufficient condition for infeasi- 
bility. In addition, a proof was provided that the value of the minimal 
feasible flow could be obtained by changing the sign on the arc length 
of the dual and finding the longest route from the destination back 
to the origin. 

During 1970, M. F. Sakarovitch was also working on these same 
problems. His research (1970) was motivated by some conversations 
with de Ghellinck. He also proved that a cycle of negative length was 
a necessary and sufficient condition for infeasibility. In addition, 
he proved that an optimal flow in a primal arc of a maximal flow net- 
work is equal to the difference between the dual shortest distances 
from the origin to the two nodes at opposite ends of the dual arcs 
intersecting the primal arc. This important property will play a 
major role in the algorithm to be presented below for non-planar 
directed networks. 

All of the work just summarized is restricted to those networks 
which are source-sink planar. The main purpose of this paper is to 
extend these ideas to general flow networks. 

In the sections to follow, we will first review the ,f state of 
the art ,f for source-sink planar directed network flow problems. We 
will then turn to the problem of non-planar networks . A pseudo dual 
will be defined and the special form of the dual shortest route problem 
will be stated. The details of the algorithm for solving this problem 
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will then be presented followed by an example, and the proofs which 
validate the algorithm. The final section presents a brief discussion 
of the usefulness of the algorithm. 
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2. MAXIMAL FLOW IN A PLANAR NETWORK 

The original flow network will be called the primal network. 

A mesh of a planar primal network is any region surrounded by nodes 
and arcs but containing neither in the plane on which the network is 
constructed. The region of the plane completely surrounding the 
primal network will be called the external mesh. 

The construction of the dual of a source-sink planar directed 
network consists of the following steps (McMasters (1970)): 

1. Connect an artificial arc between the sink and source 
of the primal and position it below the network such that it 
crosses no other arcs. The resulting network will be referred 
to as the modified primal network. 

2. Place a node in each mesh of the modified primal in- 
cluding the external mesh. Let the origin of the dual be the 
node in the mesh involving the artificial arc and the destina- 
tion be the node in the external mesh. 

3. For each arc in the primal (except the artificial arc) 
construct two oppositely directed arcs that intersects it and 
join with nodes in the meshes adjacent to it. 

4. Assign the value of the upper bound capacity of the 
primal arc as the length of the intersecting dual arc having 
the same direction that the primal arc would have if it were 
rotated 90 degrees counterclockwise. Assign to the oppositely 
directed dual arc a length equal to the negative of the lower 
bound capacity of the primal arc. 
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The convention described in step 4 is illustrated in Figure 1 

where the arc (i,j) is the primal arc having its flow X.. constrained 

iJ 

between two values L.. and M.. such that L. £ X £ M . The 

ij ij ij 

adjacent dual nodes are u and v. Arc (u,v) is assigned a length 

equal to M while arc (v,u) is assigned a length equal to -L. . 

1 J ij 

We assume M ^ 0 but place no restriction on L.. except that it 
ij ij 

cannot exceed M. . . 




Figure 2a is an example of a flow network having node 1 as its source 
and node 4 as its sink. The numbers on each are represented L_„ , 

. Figure 2b shows the initial phase of the construction of the 
dual. The dual nodes are designated as A, B, C and D. Each 
dashed arc corresponds to the location of each pair of dual arcs. Fig- 
ure 2c is the completed dual with its associated arc lengths. 

Ford and Fulkerson (1962) defined the generalized value of a 
primal cut set when lower bounds are positive by equation (1) where 
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the set S contains the primal source node and the set S contains 
the primal sink node. 



V(S,S) = l M - l_ L (1) 

ieS J ieS J 
jeS jeS 

They then extended the min-cut max-flow theorem to include 



this definition of the value of a cut set; namely, if a feasible flow 
can be found for a network with nonnegative lower bounds on arc flows 
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then the maximal possible flow through the network is equal to the 
minimal generalized value of all possible cut sets of the (S,S) type. 
Later Fulkerson (1962) dropped the restriction that the lower bound 
be nonnegative. 

The relationship between the convention of Figure 1 and the 
generalized value of a primal cut set is described by the following 
theorem: 

Theorem 1 : Any dual chain without cycle* directed fa om 

the dual *ource to the dual *ink ha* a length equal to 
the generalized value ofi the primal cut *et inter* ected 
by the arc* ofa the chain . 

Proof: Each path without circuits through the dual of a source- 

sink planar undirected graph intersects a cut set of the primal network 
such that the primal source and sink are separated by the cut and 
there is a separate path corresponding to each cut (Ford and Fulkerson 
(1956), Whitney (1933)). This is still true if the primal and dual 
arcs are replaced by directed arcs. Further, because the dual graph 
under the Figure 1 convention has two oppositely directed arcs between 
every pair of nodes in adjacent primal meshes it is also true that a 
separate chain exists directed from the dual source to the dual sink 
corresponding to each primal cut. 

Suppose we select an arbitrary cut set of the primal network 
such that the source is in a set S and the sink is in set S. Any 
such cut set can be illustrated as shown in Figure 3 for planar graphs. 
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Notice that this set contains both arcs directed from S to S and 
from S to S. Let the chain directed from the dual source to the 
dual sink which intersects the cut set arcs be a-b-c-d-e-f. 




Using the convention of Figure 1 we assign a length of -L 
to dual arc (a,b) which crosses primal arc t, a length of M g to 
arc (b,c) which crosses arc s, a length of to arc (c,d), 

a length of -L to arc (d,e), and a length of to arc (e,f). 

The chain directed from a to f would therefore have a length of 

M+M+M-L -L 
s r p t q 

which we see immediately corresponds to the form of the right-hand 
side of (1) for this arbitrary cut set (the arcs p, r and s have 
ieS and jeS while arcs q and t have ieS and jeS). Because 
our cut was arbitrary we realize that if we analyze any cut set of the 
(S,S) type in this manner we will get the length for its dual chain 
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being described by this general form. Therefore, each chain without 
cycles from the dual source to the dual sink will have a length equal 
to the generalized value of the primal cut set it intersects. 

As a consequence of Theorem 1 and the generalization of the 
min-cut max-flow theorem we get the following corollary: 

Cofiolla/iy 1: The value ofa maximal falow In a keoAible 

cLOiected 4ouAce-4ink planaA. network li equal to the 
length o i the AhoAte&t fioute through the dual network. 

Recalling Figure 2a we know that the value of the primal cut 
set consisting of arcs (1,2), (2,3), and (3,4) is 

M 12 + M 34 " L 23 = 5 + 3 " 1 = 7 * (2) 

The route through the dual network which corresponds to this cut set 
is the chain [ (A,C) , (C ,B) , (B ,D) ] which has a length given exactly 
by equation (2) . This chain happens to be a minimum route of Figure 
2c so that the maximum possible flow through the primal network is, 
in fact, 7 units. As we will show below, the fact that we found a 
shortest route of finite length means that the primal problem is feasible. 

The extension of the min-cut max-flow theorem involving the use 
of the generalized value of a cut set as given by (1) is contingent 
on the existence of a feasible flow. Feasibility requires that a set 
of arc flows can be found which satisfy the flow conservation equations 
at each node as well as the upper and lower bounds on individual arc 
flows. When the lower bound is nonpositive then we know that a flow 
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problem will always be feasible since =0 is a feasible solution 

under the M_ ^ 0 assumption. If, however, the lower bounds on some 
arcs are positive then we cannot claim such problems will be always 
feasible. An example of an infeasible problem can be obtained by 
changing the lower bound on arc (1,3) in Figure 2a to a value of 
3. 

The answer to the feasibility question is provided by Theorem 2. 



Theon.cn 2: A necet>Aan.y and buHicient condition ion. an 

IniwAiblc filou) problem aAAociated with a directed 
6 o voice- ^ ink planaA nctwon.fi i6 that a cycle o& negative 
length exJj>tt in the dual network. 

The proof of Theorem 2 depends upon the circulation theorem of 
Hoffman (1960) which is addressed to networks having a sourceless- 
sinkless structure. Such a network is easily created from any flow 
network by connecting a directed arc from the sink to the source and 
specifying the appropriate upper and lower bounds on its arc flow. 

The circulation theorem states that a network flow problem will 
be feasible if and only if 



l 

ieY 

jeY 



M. . ;> 



l L -- 

ieY 

jeY 



(3) 



for all cut sets (Y,Y) of the network after it has been converted 



to the sourceless-sinkless form. 



Suppose we define the value of the cut set (Y,Y) 



as 



V(Y ,Y) = I M ±j - l_ L. 



ieY 

jeY 



ieY 

jeY 



ij 



( 4 ) 



If we do this then the feasibility condition given by (3) is equiva- 
lent to V(Y,Y) ^ 0 and infeasibility corresponds to some V(Y,Y) < 0. 

An interesting correspondence between dual cycles and V(Y,Y) 
is described by Theorem 3. We will prove this theorem before proving 
Theorem 2 since it will be useful in the proof of the latter. 



Theorem 3 : Tf ivie Ia a Aepa/iate dual cycle InteAAecttng 

oaca o{± each pnlmal cut Act (Y,Y). TklA cycle haA a 
length equal to V(Y,Y). 



Proof: Any circuit in the dual corresponds to a cut set (Y,Y) 

of the primal (Whitney (1933)). Because of our convention of having 
two oppositely directed dual arcs intersect each primal arc we know 
that two cycles exist (oppositely directed) in the dual corresponding 
to each cut set of the (Y,Y) type in the primal. Suppose the set 
Y is circumscribed by these dual cycles. The convention of Figure 1 
will then result in the cycle directed clockwise having a length 
corresponding to V(Y,Y). If Y is the set circumscribed by the 
cycles then the counterclockwise cycle has a length V(Y,Y). 

Figure 4 is an example of an arbitrary cut set of the (Y,Y) 
type with its associated dual cycles. Primal arcs e, f, h and j 
are directed from Y to Y; arcs g, i, k and l are directed 
from Y to Y. The convention of Figure 1 will result in the cycle 
p-q-r-s-t-u-v-w-p having a length of 
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M + M.+R+M - L - L - L, - L 
g i Tc e f h j 



which is the right hand side of (4) . 




Figure 4. 



We will now prove Theorem 2. Suppose that the primal problem 
is infeasible. From the circulation theorem we know that there will 
be some cut set (Y,Y) of the sourceless-sinkless primal which has 



l M. . < 

ieY ieY 1J 






jeY 



We can rewrite (5) in the following form: 



(5) 



y m. . - y l. . < o. (6) 

ieY 1J ieY 1J 
j eY j eY 

From (4) we know that the left-hand side of (6) is the value of the 
cut set (Y,Y). It is also the length of a cycle in the dual because 
of Theorem 3. Therefore, an infeasible primal problem results in a 
cycle of negative length in the dual. 
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To show sufficiency, we consider a cycle of negative length in 
the dual. Let its length be represented by D. Because of Theorem 3, 

D = l M - l_ L = V(Y,Y), 
ieY J ieY J 
jeY jcY 

and D being negative means that the value of some cut set (Y,Y) in 
the primal is negative. For a cut set (Y,Y) to have a negative value 
(5) must be satisfied and the primal will be infeasible. 
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3. MINIMAL FLOW IN A PLANAR NETWORK 



It is sometimes useful to also know the minimal feasible flow 
for a network having positive lower bounds on arc flows. The following 
theorem provides an easy way of determining such a flow. 

Theorem 4- Ifi the anc lengths ofi the dual ofa a feasible 
pnimal filotA) netwonk have theln signs nevensed then the 
longest chain directed finom the dual destination back 
to the dual bounce Mill have a length equal to the value 
oft the minimum feasible filoM thnough the netMonk. 



The proof of this theorem is based on the min-flow equivalent 
of the max-flow min-cut theorem (Ford and Fulkerson (1962)); that is, 



min Q = max £ L , “ l_ M ,- , > ' (7) 

icS J ieS J 
jeS jeS 

where Q represents the total flow leaving the network. 

The combination of our convention of allocating dual arc lengths 
and the sign reversal just described will cause each chain from the 
destination back to the origin in the dual network to have a length 
equal to 



ieS 

jeS 



ij 



- y m. . 

ieS 1J 
jeS 



The proof of this property parallels that of Theorem 1. Therefore, 
we get the right-hand side of (7) when we find the longest chain back 
through the dual. 
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4. OPTIMAL ARC FLOWS IN PLANAR NETWORKS 



The final aspect of the source-sink planar network analysis is 
the determination of the optimal arcs flows of the primal network. 

These flows are available from the solution to the dual network shortest 
route or longest route problem depending on whether we want maximal or 
minimal flow through the network. 

Consider the maximal flow problem. Let A be the set of all 
dual arcs associated with a given primal flow network. Then the dual 
network shortest route problem can be stated as follows: 



minimize 



^ A ij 




subject to 



N 

l 

j=l 



[U.. 

Ji 




(i = 1), 

(i = 2,3,... ,N-1) , 
(i = N), 



( 8 ) 



and IK. :> 0 ((i, j)eA) . 

If we think of U.. and Z.. as representing respectively a 
ij ij 

flow and a unit flow cost through a dual arc (i,j) then we realize 
that (8) corresponds to a minimal cost network flow problem for a 
total flow out of one unit. The optimal solution will yield a chain 
flow of one unit over the cheapest route from the origin to the desti- 
nation of the dual network. This route obviously corresponds to the 
shortest route through the dual network. The dual linear programming 
problem associated with (8) can be written as 
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maximize V__ - V- , 

N 1 

(9) 

subject to V_. - ((i,j)eA). 

If we apply any of the current algorithms (see Dreyfus (1969)) 

for determining the shortest routes from node 1 to every other node 

we will obtain a number f^ associated with the length of the shortest 

route from node 1 to node i. If we let V. = f. then we will have 

i i 

an optimal solution to (9) . 

The relationship between the V values and the optimal arc 
flows in the original primal network is described by Theorem 5 (Sakaro- 
vitch (1970)). 



Thzosim 5: Conside/i a phjjnaJL aAc (a,g) ctiAzcXzd fisiom 

node, a to node, g u)kich ka6 a fi&ou) X . LeX tine, 

ap 

l nteAA ecllng dual, a/iei, be (i,j) and (j,i) with the 

length oh (i,j) being M . and the. length oh (j,i) 

ap 

being -L . I jj X t>alh>h- iet> 



x . = v. - v. 

ag j i 



( 10 ) 



ho A eveay asic (a,g) then an optimal h^lble solution 
haA been hound to the pnlmal maximal blow problem. 



Proof: Because V. - V. is a feasible solution to (9) we 

J i 

know that 



a£ 



V. - Vi 



£ M 



a e 



under our convention for assigning dual arc lengths. Thus 
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SX„SM 



ag ag 



ag* 



and therefore X is a feasible flow through arc (ct,g). 

To check for flow conservation at the primal nodes we first 

introduce a new primal arc (n,l) connecting the sink to the source 

and set its upper and lower, bounds as well as the arc flow from the 

sink to the source at the value of V — V . We then have a source— 

N 1 

less-sinkless structure and the flow conservation equations can be 
written as 

l [ X «6 - V] ' ° (11 > 



for ot “ 1,2,..*, n • 

If we next add dual arcs to intersect the primal (n,l) 
arc we will create a dual network which has a cycle circumscribing 
every primal node. If we assign arc lengths to these arcs according 
to our convention, we will be adding two more constraints to (9) . 
However, the previous optimal solution to (9) is also optimal for 
this modification. 

If we sum the V. - V. differences over all arcs of any of 
3 1 

these dual cycles we will get zero since each appears twice in 

the summation, once with a plus sign and once with a minus sign. 

Therefore, X „ = V. - V. results in the sum of the flows entering 
ag j l 

and leaving a primal node being equal to zero and (11) is satisfied. 

Finally, the optimal solution to (9) will result in V_. - 
= 1 . . for all arcs on the shortest route between nodes 1 and N 
of the dual and equation (10) will result in the flows in the arcs 
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of the primal minimum cut set (S,S) being M if aeS and 8eS 

ap 

and if aeS and geS. No flow augmenting paths will therefore 

exist between the primal source and sink and hence the existing flows 
are optimal. This completes the proof of Theorem 5. 

The reader will recall from Theorem 4 that we also can get 
the minimal flow if a problem is feasible by changing the signs on the 
arc lengths of the dual and finding the longest route from the sink 
back to the source. The associated longest route linear programming 
problem can be stated as 



maximize 7 7. I . . U.., 
L L k ij ij 



N 



+1 



subject to £ [thj - U ^ 0 



j=l 



-1 



(i = N), 

(i = 2,3,... ,N-1) , 

(i = 1) , 



( 12 ) 



and U . . ^ 0 . 
ij 

The corresponding dual linear programming problem is 



minimize - W^, 

subject to W. - W. ^ ((i,j)eA). (13) 

l j ij 

If we determine the longest route from each node to the 
origin, we will have an optimal feasible solution to (13). The value 
of is minimal feasible flow through the network (a negative 

value means flow from the sink back to the source) . The optimal 
feasible solution to the original primal minimal flow problem can then 
be obtained from the results of Theorem 6. The proof of Theorem 6 
parallels that of Theorem 5. 
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Theorem 6: Con6ideA a pfvunal oac. (a, 3) cUAected ^nom 

node a to node 3 which ha6 a ilow X . Let the 

Otp 

inteu eating dual aae 5 be (i,j) and (j,i) with the 
length o{, (i,j) being -M^ and the length o£ (j,i) 

being . I £ x Aatis&let 

a3 j i 

fioa eveny one (a, 3) then an optimal 6 ecuible solution 
hca been fiound to the primal minimal 6 low problem. 
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5 . AN EXAMPLE 



The results of Theorems 1 through 6 will be illustrated by 

the following example. Figure 5a is a network having node 1 as its 

source and node 2 as its sink. The numbers on the arcs represent the 

lower and upper bounds on arc flow. The dual construction is shown 

by the nodes A, B and C and the dashed arcs connecting them. 

The resulting dual network for the maximal flow problem is shown in 

Figure 5b. The maximal flow through the primal is 4 corresponding 

to V - V . The dual shortest route is the chain of arcs in Figure 
L* A 

5b having the hash marks. The primal minimum cut set is therefore 

{( 1 , 2 ), ( 2 , 3 )}. 




Figure 5. 



The flow through arc (1,2) of the primal is 6 as given by 

the difference V - V since the dual arcs (B,C) intersect it. 

C B 

The flow through arcs (2,3) and (3,1) are both 2 corresponding 

to the difference V. - V_ . 

A B 
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If is changed to a value of 4 then the top (A,B) 

arc of Figure 5b has a length of -4. A cycle of negative length 
will exist consisting of the top (A,B) arc and the bottom (B,A) 
arc. The problem is infeasible according to Theorem 2. Inspection 
of Figure 5a under this change shows that an infeasibility condition 
has been created by arc (3,1) having a lower bound exceeding the 
capacity of arc (2,3). 

Figure 6 is the dual network for the minimal flow problem . 



I 




The minimum feasible flow from the source to the sink of the primal 
is = 1. The arc flows creating this are “ w c ~ W B = ^ 



and X 23 = X 31 = W A - W B 



= 3. 
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6. MAXIMAL FLOW IN NON -PLANAR NETWORKS 

We would like to make use of the ideas developed in the preceding 
section for all networks and not just those which are source-sink planar. 
The purpose of this section is to present a way of approaching such 
problems . 

Suppose we construct a two-dimensional representation of a non- 
planar network such that all arcs are straight lines. If we do this, 
then no two arcs will intersect more than once at points which are not 
nodes. If we designated these extra arc intersection points as pseudo 
nodes we will have created a network which is source-sink planar and 
we can construct its dual network using the procedure and convention 
of section 2. 

Unfortunately, we have no guarantee that the shortest route 
through this dual network will give the maximal feasible flow and that 
the associated optimal feasible arc flows can be determined from 
Theorem 5. Consider Figure 7. Suppose arc (a,£) and (y>6) intersect 
at pseudo node <j> . 




Figure 7. 
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If we hope to use the results of Theorem 5 to obtain the 

optimal feasible flows in these arcs, we must require that the 

values V , V., V , and V in (9) for the associated dual nodes 
l J k m 

satisfy 




Condition (14) is needed because an optimal feasible solution to the 
original primal maximal flow problem must have 

- V 

\* ■ V 

since arcs (a, 3) and (y,6) do not intersect in the original primal 
network. Therefore, for Theorem 5 to be useful, condition (14) must 
be added to the constraints of (9) . 

We note in passing that (14) can be consolidated into only one 
equation which is 

V. + V. - V - V = 0. (15) 

j k i m 

Let the node group surrounding a primal pseudo node cf) be 

represented by (i,j,k,m) . Let $ be the set of all pseudo nodes 

<P 

in a given primal network. The incorporation of (15) into (9) then 
results in the following problem: 
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minimize V XT - V_ 

N 1 

subject to V_. - £ £ _ ((i»j)eA), (16) 

v , + \ - V. - v = 0 ( (i» j »k,ra) ; 

j k. i m 9 

Obviously, the optimal solution to (16) cannot always be deter- 
mined using an algorithm for finding the shortest route from node 1 
to each node i. The following algorithm is, therefore, suggested as 
a way of solving (16) . 

Algorithm: 

1. Compute the shortest route lengths from the dual origin to all 
nodes. If, in the process, a cycle of negative length is detected, 
terminate the algorithm. The primal flow problem is infeasible. 

If no cycle of negative length exists then compute the value of 

V. + V. - V. - V for (i,j,k,m). associated with each pseudo 
jkim J <j> 

node in $ and go to step 2. 

2. Examine the set of V. + V, - V. - V values obtained in step 1. 

J k l m 

(a) If V.+V, - V. - V = 0 for all pseudo nodes in $ then 

J k l m 

terminate. The optimal solution to the primal network 
problem can be obtained by applying Theorem 5 to the planar 
representation of that network. 

(b) If V. + V. - V. - V ^0 for any set (i,j,k,m) then 

jkim 

compute a set of numbers Y_^ (i = 1,2,...,N) where 

y n = v n and 

Y . = max (Y . -£ . . ) 

1 j J 1J 

for i = 1,2,...,N-1. Go to step 3. 
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3. For each set (i,j,k,m} for which V. + V, - V,, - V t 0 

<p j k i m 

compute the following numbers: 



‘i(*> ■ min “i j - v J - Y i>W* 

t 2 4) - max{-J jl ,Y j -V.,Y in -V k }. 

If £-^((})) < ^ 2 ^) f° r one or more sets then go to step 4; 
otherwise go to step 5. 



4- Determine that set (i,j,k,mK corresponding to 

<P 

A A 

£ 2 (<f>) - = max{«, 2 ((},) - ^( 4 ,) j (<J>) < i 2 ((|.)} 

and change z . . > Z . . , Z 1 and Z , of the arcs asso- 
i j ji km mk 

ciated with that set to 





= ^(<* 0 , 








Return to step 1. 



5. 



Select any set (i.j.k.m)^ for which £,(<}>) 

9 1 

values have been computed and change Z . . , 



Z , of the arcs associated with that set to 
mk 



and 




* 2 (<i>) 

Z, and 
km 



Z\ . = Z! = Z. (4>) 
ij km 1 Y 



Z \ . = Z\ = -JL (<J>) 
ji- mk 1 r 



Return to step 1. 
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7 . AN EXAMPLE 

To illustrate the algorithm we will solve for max Q in the 
example shown in Figure 8. This network, while being planar, is not 
source-sink planar. Figure 8 is , in fact, the two dimensional repre- 
sentation. We see that all arcs are straight lines and one pseudo 

node has been created. The numbers on the arcs are the M. . values: 

iJ 

L . . = 0 for all arcs . 
ij 

The dual network is shown in Figure 9. The and Y_^ values 

are shown next to their respective nodes. We see immediately that 

V._ + V - V~ V ^ 0 so the Y. values are needed. Next we obtain 
f e d g r 

= 5 and = Therefore, we change and i to 

5, and to -5 and recompute the shortest route. 

The new set of values are shown in Figure 10. We now have 

+ V - V. - V =0 so we are done. The maximal possible flow 
f e d g 

through the primal network is 7 as indicated by the difference 

- V^. The optimal solution to the maximal flow problem associated 
with the primal network is provided by Figure 11. The numbers on the 
arcs represent X„/M_. The minimum cut set is { (4 , 5) , (4 , 2) , (1 , 2) , 
(3,2), (5, 3)}. 



I 
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V b =9 Q V c =IO 




V h =IO 

V 10 



Figure 9. 
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0 




Figure 10. 




Figure 11. 
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8. PROOF OF THE ALGORITHM 



When V.+V. - V. - V ^ 0 for some set (i, 1 ,k,m) . we 

j k i m * J * 9 <t> 

must seek new values, V!, V ! , V. ! , V f , which will satisfy the 

r j k m J 

constraints of (16) . In an attempt to avoid reducing the value of 

and, therefore, the value of max Q we first compute the 

Y. values. The differences Y. _ - Y . (i = 1,2,...,N-1) are the 
1 N i 

lengths of the shortest routes from each of the nodes to node N. 
Now 



V N ' 



V. £ 
1 



y n ” 



(17) 



since will be less than or equal to the length of any 

chain directed from node i to node N (equality holds only if 
node i is on the shortest route from node 1 to node N) . Because 
Y n = V , (17) reduces to 

Y. £ V. . 

1 1 

Therefore, if we select satisfying the constraints of (16) such 

that 

Y. <: V! £ V. (18) 

x 1 i 



then we are assured that will not change. 

Combining (18) with the inequality constraints of (16) results 
in the following ranges for Vl - and - V^: 



max{-£ . . >Y . -V } £ V! - V! £ min{£ ,V -Y } , 
Ji J i 3 1 13 31 



( 19 ) 



"“‘■‘t-Vk-w s v ; - v k s minU km’ w • 
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Let us assume the ranges in (19) have some values in common. 

We next combine the inequalities of (19) with the equality constraints 

of (16) . The result is (20) . Notice that Z . and JL have been 

mk. km 

omitted. This is due to the fact that = l and 

ij m£ j l mk 

for every (i,j,k,m) as a consequence of the dual construction. 

<P 



max{-£. .,Y.-V.,Y -V. } £ V| - V! = V 1 - V/ £ minU . . ,V . -Y . ,V -Y }, 
ji j i m k j i m k ij j i* m k 



We realize that (20) reduces to the following form: 



V*> £V j - V i- V m - V k St l ( * ) > 

when we use the definitions of J£^((f>) and from step 2 of the 

algorithm. 

Step 4 of the algorithm handles the problems where the ranges 
of (19) have values in common for all <f> which violate (15) . It 
changes the arc lengths of arcs (i,j) and (k,m) to £,^(<f>) and the 
arc lengths of (j,i) and (m,k) to -Jl^Of)). This forces the condi- 
tions 



while satisfying 
(k,m) , will have 



V 1 - V 1 = V - V 
N INI 



V - V, = V. - V. = JL (<()) , 
m k i l 1 Y 



V. - V, = V. - V , 
i k j m* 



(18). Note that one node of each pair, (i,j) and 



= V. after the new computations. Finally, 
so that the value of maximal flow has not been 



( 20 ) 



affected by the new computations. 
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If the ranges of (19) have no values in common for some set 

(i,j,k,m) then &(<}>) < JL (4>) and the current V„ - V, value 
<f> 1 l N 1 

represents an infeasible value for max Q. Because must be 

reduced by the amount 

= max{ & 2 (4> ) - ^(ij)) | (4>) < £, 2 (<(>)} 

<f> 

before all sets (i,j,k,J.)^ can meet the constraints of (16), step 3 
adjusts the appropriate arc lengths of the set or sets corresponding 

A A A 

to to £^(<t>) and -£^(<j>) to provide this reduction in 

immediately. No further use of step 3 is necessary after this reduc- 
tion. It is possible, however, that adjustments of the step 4 variety 
will still be needed before an optimal solution is reached. 

Changing the arc lengths of the arcs (i,j), (k,m) and (j,i), 

A A 

(m,k) to and -£^(<J>) values or to and -£^(<j>) values, 

respectively, results in the optimal solution to the primal having a 

partial flow of Z^((p) or in the intersected primal arc (y , 6 ) . 

The arcs whose lengths have been changed either remain as members 

of the tree of shortest routes from node 1 to all nodes (the early 

tree) or become members in the next round of calculations. Only 

when all are members of an early tree will the condition V. - V. = 

J 1 

V - V, be satisfied if it was not satisfied initially, 
m k 

The forcing of the reverse arc of a pair to a value of -&^(<1>) 
or does not create a primal infeasibility (that is, a cycle of 

negative length in the dual) for the following reasons. Suppose arc 
(i,j) was originally on the early tree. The condition L ^ £ M ^ 
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means that the minimum possible value of Z.. will be -M r . This 

Ji Y<5 

will occur only if in the primal problem. If arc (i,j) 

was originally in the early tree then the shortest route from node i 
to node j is over arc (i,j). Any other route in the dual network 
will have a length at least as great as S, = . Thus, even 

specifying = _ ^y(S does not create a negative cycle with any of 

the routes. The worst that happens is that a cycle of zero length 
exists consisting of the arcs (i,j) and (j,i). If arc (i,j) is 
now shortened to a value of or Jl^(<j>) then it remains on the 

tree and & will have a value or -£^(<j>) which is larger 

than -M^ so no negative cycle is created. 

Next, suppose arc (i,j) was not originally in the early tree. 
There must, therefore, be some other route between nodes i and j 
which is shorter (it might be arc (j,i)). Reducing the length of 
(i,j) to £ (<j)) or &^(<f>) causes it to become an alternate member 
of the early tree. Changing the length of arc (j,i) to a length of 
or -£,(<{>) creates a cycle of zero length with arc (i,j). 

If the original shortest route between nodes i and j was not the 
arc (j,i) but rather some other route directed from node i to 
node j then its length is at least Jt^(<j>) or & (<J>) or arc (i,j) 
would not be the alternative shortest route after the new V_^ computa- 
tions. Therefore, at worst a cycle of zero length would be created by 

A 

setting & at a value of -£^(<J>) or -£^(<f>) . If, instead, the 
original shortest route was some route directed from node j to node 
i then a change in has no effect since it is not a member of a 

cycle containing the original shortest route. Reduction of to 

the value of or £^(<j>) will at worst create a cycle of zero 

length when combined with this route. 
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9. COMMENTS ABOUT THE ALGORITHM 

The intent of the algorithm of section 6 is to provide a means 
of solving dual networks for primal network optimal flows without the 
worry of whether the network is source-sink planar. All that is required 
is that the primal network be drawn such that no arc intersects any other 
arc more than once. While the algorithm is designed only to solve 
maximal flow problems, the necessary modifications for solving the 
minimal flow problem are minor and obvious. 

This algorithm may not be as efficient an approach as the maximal 
flow algorithm of Ford and Fulkerson even if no pseudo nodes are needed 
because the dual graph must be constructed. The existence of pseudo 
nodes could reduce efficiency further since several iterations of the 
algorithm may be required before an optimal flow is found or the primal 
network is found to be infeasible. Careful construction of the two- 
dimensional representation of the primal can, however, reduce pseudo 
nodes to a minimum. 

The algorithm is quite useful and efficient for sensitivity 
analyses and parametric studies of network flow as a function of primal 
arc capacities and for problems involving capacity expansion and reduc- 
tion (see, for example, Doulliez and Rao (1971) and McMasters and Mustin 
(1970)). In such problems, the dual network needs to be constructed 
only once. The changes in primal arc capacities then appear as arc 
length changes in the dual and a new shortest route can be quickly 
determined when necessary for each change. 
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